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INTRODUCTION 


Welcome. You are about to step into a field that, ten years ago, was open 
only to users of large computer systems costing tens or hundreds of thousands 
of dollars, and even five years ago, required a several-thousand-dollar calcu- 
lator that occupied the better part of a desktop. Today, the HP-25 puts pro- 
gramming into the hands of the individual. It is hoped that this book will 
allow you to realize some of the potential of this calculating instrument. 


These HP-25 Applications Programs have been drawn from the varied fields 
of mathematics, statistics, finance, surveying, navigation, and games. They 
have been arranged in eight chapters which follow roughly the above classifi- 
cation. Each program is furnished with a full explanation which includes a 
description of the problem, any pertinent equations, a list of keystrokes to be 
entered into program memory, a set of instructions for running the program, 
and an example or two, with solutions. To use the programs does not require 
any proficiency in programming, but some familiarity with the HP-25 Owner’s 
Handbook is assumed. 


For users who want to enhance their understanding of programming prin- 
ciples and techniques, a number of programs are provided to help in this 
respect. The first program in each chapter contains, in addition to the usual 
explanations, a more detailed description of the problem, a commented list of 
the program keystrokes with a step-by-step tracing of the contents of the 
stack registers, and a list of the keystrokes required to solve the example pro- 
blem. Whenever an interesting programming technique is used in one of these 
programs, it is described in a short section headed “Programming Remarks”, 
which, if present, will immediately precede the list of program keystrokes. 


Thus, whether your interest lies in solving a particular problem in a specific 
area, or in learning more about the programming power of your calculator, we 
hope that this book will help you get the most from your HP-25, 





| | 
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4 | 
A WORD ABOUT PROGRAM USAGE 


Various kinds of information are provided to explain the use of each program. 
Besides a short description of the problem, a list of applicable equations, and 
an example problem with solution, there are two forms that deserve some 
explanation: the Program form and the User Instructions form. 
Two different Program forms are provided, one of which is just a simplified | 
version of the other. The detailed form is used for a total of eight programs, 
one per chapter, with the simpler form serving for the rest. A section of a 
detailed form, taken from the Plotting/Graphing program in Chapter 1, is 
shown below: 











Use polar-to-rectangular for 





Vy = v cos @ = horiz. vel, 








vy = vsin 6 = vert. vel. 








] Initialize: t=O 





Start of loop 





Next time interval: 
tet+At 











The rightmost column, headed REGISTERS, explains what variables are 
stored in storage registers Ro through R7. The rest of the form is divided into 
eight columns. The first two columns describe the appearance of the display 
as the program is being keyed in: LINE shows the step number for the current 
instruction and CODE denotes the numeric keycodes corresponding to the 
keystrokes in the next column, KEY ENTRY. The entries in this column are 
the keys that must be pressed to enter the program into program memory. 
The key is denoted in this column as[@]; all other key designations 
are identical to those appearing on the HP-25. 

The next four columns, X, Y, Z, and T, trace the contents of the stack 
registers as they would change during execution of the program in RUN mode. 
Each entry under X, Y, Z, or T gives the contents of the respective register 
after the instruction on that line has been executed. The COMMENTS column 
contains additional step-by-step explanation of the program’s calculations. 

These last columns, X, Y, Z, T, and COMMENTS, are provided to help the 
interested user acquire a detailed, in-depth understanding of a particular pro- 
gram, or of programming techniques in general. 

The simplified Program forms contain the same information as the detailed 
forms except for the omission of columns X, Y, Z, T, and COMMENTS. 

The User Instructions form is the user’s guide to operating the program to 
solve his own particular problem. This form, which is composed of five 
columns, is illustrated below for the same program from Chapter 1, Plotting/ 
Graphing. 

















Store gravitational constant 





Input angle and initial speed 





Perform steps 5 and 6 any num- 





ber of times: Display time and 





horizontal distance 














Reading from left to right, the STEP column gives the instruction step 
number. The INSTRUCTIONS column gives instructions and comments con- 
cerning the operations to be performed. Steps are executed in sequential order 
except where the INSTRUCTIONS column directs otherwise. 

Normally, the first instruction is “Key in program”, which means to store 
the keystrokes of the program in program memory (switch to PRGM mode, 
press [f] , key in the program, then switch back to RUN mode). 

Repeated processes, used in most cases for a long string of input/output 
data, are outlined with a bold border, as in steps 5 and 6 above. In this case, 
the steps are repeated in order to generate a number of (x,y) pairs for a graph. 

The INPUT DATA/UNITS column specifies the input data to be supplied, 
and the units of data if applicable. The KEYS column specifies the keys to be 
pressed. [4] is used for the key, and all other key designations are 
identical to those appearing on the HP-25. Ignore any blank positions in the 
KEYS column. 


Some programs are complex enough that users have to press additional 
keys to generate some results. Those keys are also shown in the KEYS column. 

The OUTPUT DATA/UNITS column shows intermediate and final results 
that have been calculated either from the keyboard or from an executing pro- 
gram, and the units of data if applicable. Parentheses around an output 
variable, such as (t) in step 5, indicate that the result is displayed only briefly 
by a PAUSE instruction ([#] ). 





| CHAPTER 1 ALGEBRA AND NUMBER THEORY 


| PLOTTING/GRAPHING 


Most people who have labored through a ninth-grade algebra course probably 
| still respond with a shudder to the word “graph”. Evidently the tedium of 
finding y = 3x? -— 4x + 4, for integer values of x from - to +, has etched 
permanent memories in us all. Fortunately, we need not endure this tedium 
any longer. The HP-25 lends itself perfectly to this kind of repetitive 
calculation. 


The basic idea is to generate (x, y) pairs by keying into program memory the 
keystrokes required to calculate y, assuming x is given. Then the user need 
only return to the top of memory, enter a value for x, press [R/S], and see y 
displayed within seconds. The process may be repeated for as many values of 
x as desired. The programmer can take this process one step further into auto- 
mation by also having the calculator generate each new value of x, for 
example, by adding 1 to the old value, or, in general, by adding a specified 
increment Ax. A flowchart of the process is shown below. 







Calculate 
y=f(x) 








and (ii)Ax 





Halt to 
display (x, y) 


X+XtAX 


The program used here to illustrate this process takes a slightly different tack. 
We will consider the problem of plotting the trajectory of a stone which is 
hurled into the air with an initial velocity v at an angle to the horizontal of @. 
Neglecting drag due to friction with the atmosphere, the following equations 
describe the stone’s x- and y-coordinates as functions of the time t: 





x = vt cos 6 y=vtsind > et? 


a 








ee 
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where x=horizontal distance the stone has traveled 
y = height of the stone 
g = acceleration due to gravity | 
=~ 9.8 m/s? 
= 32 ft/s? | 
These equations differ slightly from the usual graphing function in that y is 
not expressed directly as a function of x, but instead both x and y are 
expressed as functions of a third variable t. The points to be plotted are still 


the ordered pairs (x, y); but now it is the time t which should be incremented 
by an amount At. | 


Notes: 


1. Any consistent set of units may be used. 





2. This is not a general plotting/graphing program; it merely illustrates the 
method by application to a specific problem. However, some study of 
the program listing and the flowchart should enable the user to adapt 
the method to his own application. 


Programming Remarks: 


1. The components of the velocity in the horizontal and vertical directions, 
Vx and vy, are computed in one step by a conversion of v and 6 to rect- 
angular coordinates ( [f] [+R] ). The values v, =v cos 6 and vy = vsin 0 
are returned to the X- and Y-registers, respectively. 


2. A pause ([1] [PAUSE] ) is used in this program in a very typical manner, 
to display briefly the output variable t, whose values are simple (0.25, 
0.50, 0.75, etc.) and do not need to be written down. 
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REGISTERS 


Ro At 
Use polar-to-rectangular for 
Vx =v cos 6 = horiz. vel. 

In, 2 
vy =v sin 6 = vert. vel. 
Initialize: t= 0 Ro Vx 





















































el< Ls Ls |= 


< 
x 

















Pause to display t 








Halt and dispiay x 








id display y 
back for next t 





















































TTL all 
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INSTRUCTIONS pa ouirs KEYS sree 


Key in program 





Store time interval 





Store gravitational constant 





Input angle and initial speed 





Perform steps 5 and 6 any num- 





ber of times: Display time and 





horizontal distance 





Display height 


To change @ or v, go to step 4. 





To change At or g, go to 
a 








appropriate step, store new value, 
4 cae am 





then go to step 4. 


Example: 


Plot the trajectory of a stone cast upwards with a velocity of 20 m/s at an 
angle of 30° to the horizontal. Use intervals of % second between points 
plotted. Let g =9.8 m/s”. 


Solution: 

0.25 Egy [0] 9.8 {1] 30 [4]20 [4] {R/S ]———> 0.25 (t1) 
4.33 (xi) 

{Rs |] ——_--+-_— 2.19 (v1) 

— or rO i OS i(it) 
8.66 (x2) 

: 3.78 (y2) 

{as )-—_$£_—________—_-—___—_—_—>- 0.75 (ts) 
12.99 (x3) 


{Rs |-——-© 4.74 (y3) 
etc. 








a parabola. 





Height 
(meters) 























Horizontal distance (meters) 
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Continue until y becomes negative. 
The table of these results is shown below: 
The plot of these (x, y) values is made and the stone’s trajectory is seen to be 
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QUADRATIC EQUATION 


The roots x,, x. of ax? +bx+c=0 


-b + Vb? ~4ac 


are given by X12 = 5 
a 


If D = (b? - 4ac)/4a? 


is positive or zero, the roots are real. In these cases, better accuracy may 
sometimes be obtained by first computing the root with the larger absolute 
value: 


If 


In either case, = : 
X, a 


If D <0, the roots are complex, being 


-b, V4ac — b? 


utiv= 
2a 2a 
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cope | ENTRY 


ao 
74 [R/S 
























































= 


TAAAAANAGRGRAA 


|| x2y 
aa 
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INPUT OUTPUT 


Key in program 


initialize 
Enter coefficients and display D 
If D > 0, roots are real 
or 
If D <0, roots are complex of 
form u + iv 
For new case, go to step 3. aa 


Example: 



































Find solutions to the three equations below: 
1 x?+x-6=0 
2. 3x? 4+2x-1=0 
3. 2x? -3x+5=0 
Solutions: 
1. D=6.25 

XxX; = -3.00 

X2 = 2.00 

D=0.44 

X, =-1.00 

X2 = 0.33 

D=-1.94 

X1,X2q = 0.75 + 1.391 





Chapter 1 Algebraand Number Theory = 15 


COMPLEX ARITHMETIC, +, -, x, + 


Let a, t+ ib, and a, + ib, be two complex numbers. The arithmetic operations 
+, -, x, + are defined as follows: 


1. +, addition 


(a, + ib,) + (a2 + ib2) =(a; + a2) +(b; +b2)i 


-, subtraction 
(a, + ib) —(a, + ibz) =(a, - a.) + (b; —b2)i 


x, multiplication 


(a; + iby) x (a + ib.) =r, 12 elf, + 62) 


+, division 


(ai+ibi) on 
(a2 + ib.) 


el: -9) a, + ib, #0 


where r, e!: is the polar representation of a, + ib, andr, e!9? is the polar 


representation of az + ib2.In each case let the answer be x + iy. 


After a calculation is finished x is stored in Rg as well as the X-register and y 
is stored in R, as well as the Y-register. In this way arithmetic operations can 
be chained together. 
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DISPLAY KEY 





















































2402 | STO2 
= 








a ‘ . 
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Key in Program 





Store first complek number 








Key in next number 








For addition 





or 





subtraction 





or 





multiplication 





or 








For imaginary part 





For next calculation in chain, go 





to step 3. 





For new case, go to step 2. 





Examples: 
1. (1.2 +3.7i) - (2.6 - 1.91) =-14+5.6i 


2. =" = 0.25 +0.64i 


(3 + 4i) + (7.4 - 5.61) an 
3. [ae [3.1 + 4.6i] = 3.61 + 7.163 
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COMPLEX FUNCTIONS |z|, z?, '/z,./z 


A complex number z = a + ib has polar representation re. The formulas used 
to evaluate the given functions are as follows: 


lzi=r 


22 =17? ei20 


n= =e? 2#0 


Vi =* (Vt oi) =2 (xt iy) 


The answer is represented by x + iy. 







































































oes jt | oe |e 
OO nN a 
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Examples: 


1. {12 — Si] = 13.00 


2. (6-i)* = 35.00 - 12.00i 
1 
2+ 5i 


4. VJ3¥4i =+(2.00 + 1.00i) 


= 0.07 - 0.17i 
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DETERMINANT AND INVERSE 
OF A2x 2 MATRIX 


412 : 
be a 2 x 2 matrix. 
a22 


The determinant of A denoted by Det A or |A| is evaluated by the following 
formula: 


Det A =a Ay, — 442 42, 


Also, the program evaluates the multiplicative inverse A~' of A. The following 
formula is used: 


sein a2/Det A -a,,/Det A 
-a,,/Det A a;;/Det A 


DISPLAY 


LINE 
2400 |RCLO 
24 04 | RCL4 la ic 




















es 2401 | RCL1 1300 | GTO 00 









































2400 | RCLO 
CA Vis 
74 | R/S 
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STEP INSTRUCTIONS 


Key in program 





Store matrix 














Compute determinant 





Compute inverse 














For new case, go to step 2. 


Example: 
Find the determinant and inverse of the matrix 


[2] 


Solution: 


Det A = -20 


acta [920 0.10 
0.20 -0.15 
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NUMBER IN BASE b TO NUMBER IN BASE 10 


This program consists of two subprograms. The first changes the integer part 
of a number in base b to a number in base 10. 


x rer 2 inady. 3 n-2 i i 
lho = lnlp-y - doh =i,b" i Fip_i b +...+ib +i; 


This is evaluated in the form 


b & (b (b (in bt in-1) x in-2) + sea) + iz) + iy 


The second subprogram changes the fraction part of a number in base b to a 
number in base 10. 


Fo =fibuw tp =f! +b +485" 


Together the two programs can convert any number in base b to a number in 
base 10. Zeros must be entered in their proper place. 
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Input final digit 





For fractional part, input digit 





after decimal 














* The stack must be maintained 





at these points. 





Examples: 
1. 17773 = 10234 


2. 143.2044; = 48.4384,, 
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NUMBER IN BASE 10 TO NUMBER IN BASE b 









This program will convert any positive number in base 10, Nj, to a number 
in base b, Ny, where 2 < b < 100. The algorithm used is an iterative one 
which adds one more digit to N, at each iteration. The program pauses as 
each new Ny is computed to display successive approximations to the final 
answer. When the displayed value of Ny has reached the accuracy desired by 
the user, he should press to halt the program, then [3] to display Np. 








Notes: 





1. When the base b is such that 11 <b < 100, two display positions are 
allocated to each digit of N,. Begin partitioning to the right and to the 
left of the decimal point. For example, 41106.12 in base 16 stands for 
4B6.C. 







An error indication during execution means that the machine’s accuracy 


has been exceeded. The value of Ny is in R3. 
REGISTERS 
Rob 


Ry Nio 
R» 10 or 100 
R3N 


R41 digit 
Rs 
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38 


3 a 
a 
== 
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INPUT OUTPUT 
a 


Key in program 
Set display format 
}3 | Store base and decimal number 
fs | When number is shown with 
| | desired accuracy, press to 
halt, then 
For new case, go to step 3. 


Examples: 


1. 67.3249 = 403.05011446 
= 43.51Ei6 





























mt = 3.14159265419 = 11.00100100, 
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VECTOR CROSS PRODUCT 


If A = (a,, a2, a3) and B = (b,, bz, b3) are two three dimensional vectors 
then the cross product of A and B is denoted by A x B and is calculated as 
follows: 


AxB= as * ie pes = (a, b; -a3 b b, -a, b b.- b,) 
x ‘ 3 a a ,a a gia a 
she bs ves 2 93 3 92,43 0) 1 93,4; 02-42 Dy 


Let the solution be represented by (c,, ¢2, ¢3). 
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INPUT OUTPUT 
ce Key in program 























Compute cross-product 











For new case, go to step 2. 


Example: 
Let A=(2,5,2) 
B= (3, 3, -4). 


Solution: 
Ax B=(-26, 14, -9) 
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ANGLE BETWEEN, NORM, 
AND DOT PRODUCT OF VECTORS 


Let a= (a,, a2 ...,a,) and b= (b;, bo, ..., by) be two vectors. 


The norm of a is denoted by lal and is calculated by the following formula: 


lal=Vay? ta,?+...+a,2 
2 
lb] =Vb,? + b2? +... +b, 


The dot product of a and 6 is denoted by a + 6 and is calculated by the 
following formula: 


similarly , 


a-b=a, b, +a, by +...t+a, by 
The angle between a and b is denoted by @ and is calculated by the following 


formula: 
@=cos7? = : . 
lal + Ib| 


The angle is calculated in any angular mode. When calculated in degrees, 
decimal degrees are assumed. 























23 51 00 


















































Chapter 1 Algebra and Number Theory 29 


INPUT 
INSTRUCTIONS DATA/UNITS. 


Key in program 





Initialize 


Perform for i = 1,..., n: 
Key in aj and bj 




















Example: 


Let A=(2,5,2) 
B = (3, 3, -4) 


Solution: 


lal = 5.74 
[bl = 5.83 
a+b = 13.00 
0 = 67.16° 
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SIMULTANEOUS EQUATIONS 
IN TWO UNKNOWNS 


Let ax+by =e 
and cxtdy=f 


be a system of two equations in two unknowns. Cramer’s Rule is used to find 
the solution. 





If ad - be = O the calculator displays Error. In this case no solution or no 
unique solution exists. 


REGISTERS 
















































































Example: 
5x -3y = 12 
2xt+y=9 


Solution: 
x=3.55 
y=1.91 
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CHAPTER 2 FINANCE 


Because many of the finance programs have certain quantities in common, a 
word about these variables and the names used to refer to them may be 
helpful. 


Five main variables recur in finance problems: n, i, PMT, PV, and FV. The 
first of these, n, denotes the total number of periods. The periodic interest 
rate i must be expressed in these programs as a decimal. Thus an annual 
interest rate of 6% is expressed as 0.06, which as a monthly rate would be 
0.06/12 = 0.005. PMT refers to the amount of the periodic payment. The 
present value, PV, is the value occurring at the beginning of the first period, 
while the future value, FV, is the value at the end of the last period. 


MORTGAGE LOAN 
ACCUMULATED INTEREST/REMAINING BALANCE 


Remaining 
balance 






PV 


Remaining balance 
after payment K 





Payment 
number 





“oa. eer 


As one enters into the realm of financial calculations, one of the most striking 
revelations is how much of the repayment of a loan goes to interest. A new 
homeowner, for example, sends off his first monthly installment of $220.13 
toward repayment of a 30-year, $30,000 mortgage assumed at 8% annual 
interest. With a proud sigh and a swelling chest, the homeowner mentally 
checks $220 off the $30,000 and figures he’s well on his way. Right? Well, 
not quite. In fact, $200 of that payment will go to interest, and only $20.13 
to reducing the principal of the loan. 
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This program will allow the user to calculate the amount paid to interest, for 
one payment or over a number of payments, as well as the amount of principal 
still unpaid, i.e., the remaining balance. The user must input the following 
values: the initial amount of the loan, the periodic interest rate, and the 
periodic payment amount. He must then key in a beginning payment number, 
J, and an ending payment number, K. The program will compute the accumu- 
lated interest charge from payment J through payment K, inclusive, and the 

balance remaining after payment K. If one wishes to find the amount of 

interest paid in a single payment, he can simply set K = J. 





The program can also be used to generate a limited amortization schedule 
showing the balance remaining after successive payments. This can be done by 
leaving J = 1 and increasing K by 1 at each iteration. Outputs will be the total 
amount paid to interest over the first K payments, and the balance remaining 
after payment K. 


Equations: 





\-K 
BAL, = eee esr +p" =1 + r| 
(1+ i)-* , 


Inty_x = BAL, - BALy_; +(K-J+ 1) PMT 


where BAL, = remaining balance after payment n 


Inty_« = accumulated interest, payments J through K 


PV = initial loan amount 
PMT = periodic payment amount 


i= periodic interest rate 


Notes: 
1. The periodic interest rate i must be entered as a decimal. For example, 
for monthly payments with an annual interest rate of 9%, the periodic 


interest rate should be input as i= ae = 0.0075. 


2. The use of this program is not restricted to mortgage loans, but applies 
equally well to any loan which is being repaid with equal periodic 


payments, 








| 
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Programming Remarks: 


In many finance programs, the expressions (1 + i) and (1 + i)” are used several 
times per program. It is often simpler to calculate the quantity once and then 
store it for later use, rather than calculate it anew each time. In this program, 
the values of (1 +i)~* and (1 +i)~? are calculated once and then stored in R7, 
thus saving both program steps and execution time. The same principle, of 
course, applies to other expressions in other problems. 


REGISTERS 





Catculate BALK. 






























































BALK 





Calculate BALy_; 








R 
R 
R 
PMT s 
PMT s +PV 


o 
i 
1 
PMT 
2 
RPV 
Lets = [(1+i)7 + 
Re 
Ry i 
























































Diff = BAL - BALy-1 
K -J+ 1 gives no. PMT’s 
from J through K 




















m=K-J+1 
m PMT is $ paid, J~K 
Display Inty_« 

















Display BAL 
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INSTRUCTIONS INCU SATANINS 











Key in program 





Store the following variables: 





Periodic interest (decimal) 





Periodic payment 
4 











Initial loan amount | 


Starting payment number 


S neconedin | - (i wn ft 
~ [Rechacon apy valliblessioe™ - i nas nee eee Geel ee 


register and go to step 3. ae ay ee ‘eee 









































Example: 


A mortgage is arranged so that the first payment is made at the end of 
October, 1974 (i.e., October is payment period 1). It is a $25,000 loan at 8% 
with monthly payments of $200. What is the accumulated interest for 1974 
(periods 1—3) and for 1975 (periods 4—15) and what balance remains at the 
end of each year? Also, generate a schedule of interest paid and remaining 
balance for the first 5 years of the mortgage (periods 12, 24, 36, 48, 60). 


Solution: 


(Notice that i must be entered as a decimal, monthly rate.) 


08 [4] 12 =) BY 5) 200 EY [2} 25000 Kg [3] 1 


(4]3 EOD (5) (1) Lerem J —___—___——-» 499.33 
(interest paid in 1974) 


[Rs] —A AA __ i ——_——-» 24899 33 
(remaining balance at end of 1974) 


4 Bg (4) 15 Gey (5) (Rs) —_—_ > 1976.65 
(interest paid in 1975) 


a 2447598 
(remaining balance at end of 1975) 





ES el 
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Now, generate the amortization schedule: 


1 BE) 4) 12 BS |) ks) ——————_————__-> 1985.00 
(interest thru 15* year) 











—- > 4585.00 
(remaining balance after 15 year) 
24 Sy &) (es) > 3935.56 | 
(interest thru 2"¢ year) 
$$ 94 135.56 


(remaining balance after 2% year) 





36 By [5] ———_ > 5848 81 
(interest thru 3° year) 
25 048:81 
(remaining balance after 3°? year) 


48 Gi [5] | _—_—> 7721.67 
(interest thru 4” year) 













© 23121.67 
(remaining balance after 4** year) 
60 23 &) [es] ——-— $$? > 950.77 
(interest thru 5*” year) 
[Ris] 22550.77 
(remaining balance after 5'™ year) 
Remaining 
balance (K$) 
a8 PV = $25000.00 
i = 8% 
PMT = $200.00 
20 
15 
10 
5 
ry 


45 30 135 180 225 270 | 
Payment number 
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MORTGAGE LOAN 
PAYMENT, PRESENT VALUE, NUMBER OF PERIODS 


PV 


1 2 3 n-2] n-1 n 
Seemann ° bn. 
T PMT PMT PMT PMT 


PM 


For a loan which is being repaid with equal periodic payments, this program 
will calculate the payment amount, the present value, or the number of 
periods of the loan, given the periodic interest rate and the two other variables. 


Remember that the periodic interest rate i must be expressed as a decimal, 
e.g., 6% is represented as 0.06. 


The equations used are as follows: 


PMT = PV | —————— PV = PMT 
1-(i+i)” 1 


i t+)" 


__ In(1 -iPV/PMT) 


=. 


Ld 
oS 
oO 
Np 











1407 |fLN 
2402 | RCL2 
51 + 

39 1407 |fLN 





| 25 | 
| 26 | 
[ 27 | 
| 28 | 
| 29 | 
| 30 | 
| 31 | 
| 32 | 
38 











2402 | RCL2 








me] et | ek | oe 
On ao) a 


= 











0 1 
5 + 
3 

f 


1 
1 
2401 | RCL1 
2 | CHS 
14 03 y* 
41 
1 


24 02 | RCL2 
7 ee 








38  Chapter2 Finance 


UTPUT 


| 1 | Keyi in program 





For payment 



































For new case, go to step 2, 3, or 





4. 


Examples: 


1. What monthly payment is required to amortize a $3000 loan at 9.5% 
(.095) in 36 months? 


You are willing to pay $175 per month for 24 months on a 9.5% loan. 
How much can you borrow? 


How many months will it take to pay off a $4000 loan if your monthly 
payment is $200 and the annual interest rate is 9.5%? 


Solutions: 
Divide 0.095 by 12 to find the monthly interest rate expressed as a decimal. 
1. $96.10 


2. $3811.43 
3, 21.86 months 
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MORTGAGE LOAN 
INTEREST RATE 


PMT PMT PMT PMT PMT 


This program will calculate the interest rate on a loan with equal periodic 
payments. The user must specify the number of periods, the present value or 
initial loan amount, and the payment amount. 


The program performs an iterative solution for i using Newton’s method: 
fix) 
(ix) 
1-(1 +i)" FV 
i PMT 


ict. =i_ - 


where f(i) = 


The initial guess for i is given by 


01 
02 
03 
04 
05 
06 
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INPUT OUTPUT 
we] weer [ine | ed a 
1 


Key in program 





Store number of payments 





Key in present value and pay- 





ment amount 





PV/PMT 
Compute interest i (decimal) 


For new case, go to step 2. 











Example: 
You recently obtained a $2500 car loan for 36 months. If your monthly pay- 
ment is $86.67, what is the annual percentage rate? 


Solution: 
15.01% 
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COMPOUND AMOUNT 


FV 
1 2 3 n-2 | n-1 n 
+a Rn 
PV 


This program applies to an amount of principal that has been placed into an 
account and compounded periodically, with no further deposits. The impor- 
tant variables in this case are the number of compounding periods n, the 
periodic interest rate i, the principal or present value PV, the future value of 
the account FV, and the amount of interest accrued I. Any of these may be 
calculated from the others by these formulas: 


_in(FV/PV) 7 aa 
ip Nef m2 PV FVGtD 


FV=PV(1+i" 1=PV [(1+)" -1] 


C 
Ag|s 
Z = 
ie 

SZ 
Z 


KEY 





BP) 


| 24 06 | 

| 2404 | RCL4 | 

a ae 
1407 


| 2402 | RCL 2 
01 
51 
1407 | FUN 


R 
1 
; 























RCL 1 


= = 





























DI 
| oo | 
jot | 
| 02 | 
|_03 | 
| 04 | 
| os | 

06 
07 
z 
| 10 | 
| 44 | 
| 12 | 
13 | 
[14 | 
| 15 | 
ie 
17 
18 
19 
20 
21 
[24 | 
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Rey i in program 








To compute number of periods 





Pv 





FV 








To compute periodic interest 

















To compute principal 





i (decimal) 





Fv 








fs | To compute future value 


i {decimal) 





PV 








To compute accrued interest n 





i (decimal) 














Examples: 


1, 


Assuming an annual inflation rate of 10%, how long will it take prices to 
double? (Suggestion: let PV = 1, FV = 2) 


Find the rate of return on $1000 compounded quarterly if it amounts 
to $1500 in 5 years. 


How much will you need to invest today at 5 3/4% compounded 
quarterly to have $3000 in 5 years? 


What is the future value of $2000 invested at 5 3/4% compounded 
quarterly for 4 years (16 quarters)? 


How much interest do you receive on $1500 deposited for 10 years if 
interest at 5 1/2% is compounded annually? 





Solutions: 


1. 


2 
3: 
4 
5 


7,27 years 

.0205 quarterly = 8.19% annually 
$2255.02 (i = 0.0575/4) 
$2513.08 (i = 0.0575/4) 
$1062.22 (i= 0.055) 
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PERIODIC SAVINGS 
PAYMENT, FUTURE VALUE, NUMBER OF PERIODS 


This program calculates payment, future value, or number of time periods for 
a schedule of periodic payments into a savings account, given the interest rate 
and two of the three other variables. Remember that i must be input as a 
decimal, e.g., 6% is expressed as 0.06. 





Then n, PMT, or FV may be calculated from the following formulas: 


FVi . 
eee, + * 
7 i [Fe is | apes ee 
In (1 + i) a+! -da+i 


pV=sut [a +p) 2a +i) 


REGISTERS 


DISPLAY 


1473 | f LASTx 

2401 | RCL1 

24 02 — 1300 | GTO 00 
RCL 3 


ian 24 03 


51] + 
2300 | STOO 
51] + 















x 
RCL 3 








-_ 








1407 | f LN 





nN 
e=S 
o 
a“ 
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: 
2402 | RCL2 
01} 1 
= 
ral 
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Key in program 





To compute number of pay- 





ments 














To compute periodic payment 














To compute future value 





i (decimal) 














For new case, go to step 2, 3, 


Examples: 


i 


How long will it take to save $15,000 if you are making quarterly 
deposits of $400 at 6% annual interest? 


2. You will need $10,000 in 7 years. How large a monthly payment do 
you need to make if the annual interest rate is 6 1/2%? 

3. | How much money will a person have if he deposits $150 at the end of 
each month for a period of 3 years? He receives 6% annual interest. 

Solutions: 

1. 29.62 quarters or 7.40 years (i = .06/4) 

2. $93.82 (n = 84, i= .065/12) 

3. $5929.92 (n=36,i=.06/12) 


gE a 
| 
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DISCOUNTED CASH FLOW 
NET PRESENT VALUE, INTERNAL RATE OF RETURN 


The primary purpose of this program is to compute the net present value of a 
series of cash flows. In general, an initial investment Vo is made in some enter- 
prise which is expected to bring in periodic cash flows C,, C2, ..., Cy. Given a 
discount rate i, which must be entered as a decimal, then for each cash flow 
C;,, the program will compute the net present value at period k, NPV,. A 
negative value for NPV, indicates that the enterprise has not yet been pro- 
fitable. A positive NPV; means that the enterprise has been profitable, to the 
extent that a rate of return i on the original investment has been exceeded. 


The program may also be used iteratively to calculate an internal rate of 
return. The objective here is to find the discount rate i which will make the 
final net present value, NPV,,, equal to zero. The procedure, then, is to store 
Vo and a first guess at the rate of return i, input the cash flows C, through 
C,; and thus find NPV, If NPV,, is negative, the estimated rate of return was 
too high; if NPV,, is positive, the estimate for i was too low. Adjust the esti- 
mate for i accordingly, store the new i, and input the cash flows again. Inspect 
the new value of NPV,, to obtain a new estimate for i and repeat the process. 
The entire procedure is repeated until NPV,, is zero, or very close to it. The 
last value of i input is then regarded as the internal rate of return. 





Each figure for net present value is found by 


k C; 
NPV, =-Vo + } : 
j=1 (1+ ip 
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INPUT OUTPUT 


Store tnitial investment and 





discount rate 





Perform for k = 1,..., 9: 





Input C,and compute NPV, 


eames foe do 
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Example: | 


You have been offered an investment opportunity for $150,000 at a capital 
cost of 10% after taxes. Based on the following cash flows, will this investment | 


be profitable? 
Year Cash Flow 
1 $30,000 
2 26,300 | 
3 50,000 
4 55 ,600 
5 45,200 
Solutions: 





Remember to enter i as 0.10. 


NPV, = -$122,727.27 
NPV, = -$100,991.74 | 
NPV; = -$63,426.00 

NPV, = -$25,450.45 
NPV; = $2,615.20 


Since Cs is positive the cash flow is profitable to the extent that the cost of 
capital is 10%. 





Chapter 2 Finance 49 


CALENDAR 
DAY OF THE WEEK 
DAYS BETWEEN TWO DATES 


This program will compute the day of the week for a given date, or the num- 
ber of days between two dates, for any dates from March 1, 1700, to 
February 28, 2100. The program works by assigning the number 1 to March 
1, 1700, and a corresponding number to each succeeding day. When comput- 
ing day of the week, a O represents Sunday, 1 Monday, 2 Tuesday, etc. 


Thus for month m, day d, year y, the number N assigned to that date is 


Nim, d, y) = [365.25 g(y,m)]+ [30.6 f(m)] + D - 621049 
where 


g(y, m) = 


m+13ifm=1or2 
mt+lifm>2 


y-lifm=1or2 


yifm>2 ad =| 


[m] represents the integer function, [*][INT]. E.g., [6.34] = 
Note: 
For days from March 1, 1700, to February 28, 1800, 2 days must be added to 


the value for N calculated by the program. For days from March 1, 1800, to 
February 28, 1900, 1 day must be added. 


DISPLAY 





2401 
14 41 
13 09 








olo|l~w 
= |D)/o 





jo 
a 
je __| 


2402 | RCL2 
51 |+ R, Temporary 


06 |6 


RCL 1 
f x<y 
GTO 09 














o nip ly £ 
oa ©O|a\|N 4 














_ 
a 














1501 |g FRAC 
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a 
Key in program = i = L seems 


day 











year 





Compute N(m, d, y) 





For day of week, go to step 8 





For days between dates, store 








first N 





| Repeat steps 2 and 3 for second 





date, then 





For new case, go to step 2. 





For day of week (0 = Sunday) 











For new case, go to step 2. 


Examples: 
1. What day of the week was July 4, 1776? 


2. Find the number of days between March 27, 1948, and April 7, 1975. 


Solutions: 
1. Thursday (4). (Remember to add 2 days.) 


9872. 











(os = “sea” «wl 
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CHAPTER 3 GAMES 
MOON LANDING SIMULATOR 


Imagine for a moment the difficulties involved in landing a rocket on the 
moon with a strictly limited fuel supply. You’re coming down tail-first, free- 
falling toward a hard rock surface. You'll have to ignite your rockets to slow 
your descent; but if you burn too much too soon, you'll run out of fuel 100 
feet up, and then you'll have nothing to look forward to but cold eternal 
moon dust coming faster every second. The object, clearly, is to space your 
burns just right so that you will alight on the moon’s surface with no down- 
ward velocity. 





The game starts off with the rocket descending at a velocity of 50 feet/sec 
from a height of 500 feet. The velocity and height are shown in a combined 
display as -50.0500, the height appearing to the right of the decimal point 
and the velocity to the left, with a negative sign on the velocity to indicate 
downward motion. If a velocity is ever displayed with no fractional part, for 
example, -15., it means that you have crashed at a speed of 15 feet/sec. In 
game terms, this means that you have lost; in real-life, it signifies an even less 
favorable outcome. 


You will start the game with 120 units of fuel. You may burn as much or as 
little of your available fuel as you wish at each step of your descent; burns of 
zero are quite common. A burn of 5 units will just cancel gravity and hold 
your speed constant. Any burn over 5 will act to change your speed in an 
upward direction. You must take care, however, not to burn more fuel than 
you have; for if you do, no burn at all will take place, and you will free-fall 
to your doom! The final velocity shown will be your impact velocity (general- 
ly rather high). You may display your remaining fuel at any time by recalling 
R.. 


Equations: 
We don’t want to get too specific, because that would spoil the fun of the 


game; but rest assured that the program is solidly based on some old friends 
from Newtonian physics: 





1 = 
X= Xo t Vot + at? V=Vo + at v? =vo2 + 2ax 


where x, v, a, and t are distance, velocity, acceleration, and time. 


Notes: 


1. ‘If youccrash before running out of fuel, the crash velocity shown will be 
the velocity before the burn, rather than the impact velocity. | 


2. Use only integer values for burns. Any decimal entry will cause an error 
in the display for V.X. | 


Chapter3 Games 53 


Programming Remarks: 


An interesting feature of this program is the simultaneous display of both 
speed and altitude (V.X), as for example, -50.0500. This is accomplished by 
storing the speed V and the altitude X in their normal form (~50.00, 500.00), 
then dividing X by 10,000 (10*) before combining them. An additional sub- 
tlety involves the question of the sign of V, and whether (X/10*) is to be 
added to or subtracted from V. For example, if V=-50 and X = 500, we 
should subtract: V - (X/10*), in order to generate a display of -50.0500. But 
if V=10 and X=50, we should add: V=V + (X/10*) in order to display 
10.0050. Inspection of the program listing, lines 2 through 12, will reveal 
how a conditional branch was used to resolve the dilemma. 





Four-place display 
Form display V.X 

















Is V negative? 

Yes, branch 

No, add V and X Rg Acceler- 
ation 











V<0, add V and -X 








V.X is V + (X/10") 
Burn B has been input 
Burn > Fuel? 

















STO -2 











@)a!a]a\ n/n 
1 
a 


ill ii 
































Is X below ground? 
Yes, you've crashed 
2403 | RCL3 No, update V 
235101 New velocity; Ve V+A 
GTO 02 Display V.X 
v All fuel gone. show 
crash velocity as 
V = (V? + 29x)"/2 
where g = gravity = 5 
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INSTRUCTIONS 





DATA/UNITS 


Key in program 


2 | initialize x [soo |[ sto ][ o |[ || 500.00 
[sy initial V.X as 50,0500 | 
4 | Key in burn, compute new speed 


sraldotincs | em |Les IT IC IL. | vx 
Fe \eertonsten4 ailgoutander | as a ae oe 
6 | To see remaining fuel at any = aes eae eee 
— ; + | 
| 7 |To display speed and distance ae ee Seen oe | 
at any time (+ |[erew Jas JL___]|_ vx 
= es ee a 





















































To start a new game, go to step 2. 














Example: 


500 EE) (5) 50S SE) (7) 120 BE) I 
[F] Ceniw ] [5] 50.0500 


0 OO 55.0448 


5 -55.0393 








(note constant V when burn = 5) 


30 ————<—<—<——— a =30,0350 


0 > -35.0318 


SS 2.0000 
(remaining fuel) 

(4) [es] ——_————_____________——»  -50.0190 
(display V.X again) 

10s} ——— > 45.0143 

0  -50.0095 

(2) ——_——_—_——_—_—_————>_ 75.0000 

10 —_——_—_— > -45.0048 


25 -25.0013 


20 SS 0: 





Oops. 
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NIMB 


The game of Nimb begins with a collection of N objects, or as the calculator 
plays it, with the positive number N. Each player alternately subtracts one, 
two, or three from the total until only one is left. The player forced to take 
the last one loses. 


To begin the game, you must tell the machine how many objects to start with, 
i.e., the value of N. A reasonable number is 15. After each move the machine 
will display the remaining total. A negative sign indicates that it is the user’s 
move next, while a positive display indicates that it is the HP-25’s move. 


As the challenger you are allowed to make the first move. It is possible to win 
but of course the HP-25 is a master player: it will not let you make an error 
and win, (Not, that is, unless you cheat and take a number other than 1, 2, or 
3—a contingency so far beyond the realm of the HP-25’s naive faith in human- 
kind that the unsuspecting calculator has no way of knowing if you do or 


Ro Total 

R, Machine move 
R, + Total 

Rs 55178 

R,4 3507.1 
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OUTPUT 
DATA/UNITS 


Key in program 





Initialize 








Store total number of objects 





(usually 15) and set display 








If number in display is negative, 





key in your move Your move 





If number in display is positive, 





let HP-25 move 





Perform steps 4 and 5 until game 





is over 





At end of game, turn calculator 





upside down to read message 





For another game, go to step 3. 


Example: 
Perform the initialization with N = 15. 


User takes 3. 

3 eo «12. 

9. 
HP-25 takes 3. 


User takes 2. 

2 i 

Sr 
HP-25 takes 2. 


User takes 3. 

3s) ——— 2. 

> ele 
HP-25 takes 1. 


User takes last 1. 
1 —_—_---———__-. —_— — 55178. 


Turn calculator upside down for message (BLISS). 
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TEACH ARITHMETIC 


We at Hewlett-Packard feel that the hand-held calculator, far from threatening 
the traditional tenets of a sound mathematics education, may be used crea- 
tively to reinforce learning in such areas as arithmetic, algebra, geometry, 
trigonometry, calculus, and numerical analysis. This program, which is 
designed to be used in teaching children the four operations of elementary 
arithmetic (+, -, x, +), demonstrates some of the (largely unexplored) poten- 
tial of the HP-25 as an educational tool. 


The basic flow of the program is to pose a problem in arithmetic, check the 
answer that the user keys in against the correct answer, and then do one of 
two things: if the user’s answer was correct, the program will go on to pose 
a new problem; if the keyed-in answer was wrong, the program restates the 
original problem to give the learner a second chance. 


To run the program, the user must store a value called Max in Rg. This tells 
the program not to use any numbers as large as Max in its problems. If you 
specify a Max of 12, for example, then all the problems will deal with num- 
bers between 0 and 11. The user must then store in R, a “seed” s, a number 
between 0 and 1, which will determine the sequence of problems that will 
appear. Different seeds generate different problems, thus ensuring that the 
learning game doesn’t get boring. With the display format set to[#] [Fix] [2], 
the execution of the program will cause the first problem to be displayed as 
follows: the display will show one number to the left of the decimal place, 
and one number to the right. For example, the numbers 8 and 2 would be dis- 
played as 8.02. The user may then choose what operation to perform on the 
two numbers: he may add (8 + 2), subtract (8 - 2), multiply (8 x 2), or divide 
(8 + 2). After he keys in his answer and re-initiates program execution, the 
program will either display a new problem, if his answer was right, or display 
the same two numbers again, but this time with a negative sign in front 
(-8.02). The negative sign is an indication that the answer was incorrect, and 
does not denote a negative number. (All numbers in the problems are positive, 
though of course the results of some subtractions may be negative). If the 
problem reappears with a negative sign, the user should key in a different 
answer and try again. As soon as the correct answer is given, the program will 
go to display a new problem. 
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R, Left # 
R3 Right # 
R, Problem 












































wall |x 
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Key in program 
Store Max (0 < Max < 100) 








Store seed (O<s <1) 





Set display format 





Generate a problem 





Choose an operation and key in 





your answer: 





For addition (+) 





For subtraction (-) 





For multiplication (x) 





For division (+) 





If you were right, program will 





display new problem; go to step 
6. 








If you were wrong, program will 





|. show same problem again; go to 





step 6 again. 





Repeat steps 6—8 as many times 





as desired 





To change Max, go to step 2, 





then to step 5. 


Example: 
Let Max = 12 and the seed s = 0.725 





Solution: 
[+] [prem | [Rvs] 6.01 
(6+ 1=7) 


78) ——_———__-—__________—_—> 8.03 

(8 x 3=25) 

25 (3) (6) es) ——_—_ —————> -8.03 

(Try again: 8 x 3 = 24) 

24 GS (3) (6) (2s | ——_—_____________—_ 3.11 

(3 - 11 =-8) 

3 GS GD G8) @) es] —————_____» 9.00 

(9+0=9) 

9s |] — + 2.05 
etc, 


59 











| 61 


| CHAPTER 4 NAVIGATION 


COURSE PLANNING—GREAT CIRCLE PLOTTING 
| AND RHUMBLINE NAVIGATION 


: 4 Long voyages by sea or air are generally made to follow one of two sorts of 

routes: a rhumbline or a great circle. The rhumbline is the path of constant 
| heading between two points on the earth’s surface; it intersects all lines of | 
longitude at the same angle. It is also the course defined by the straight line 

between two points on a Mercator projection. It is a convenient course for 
| navigation because its direction does not change, and for short distances at 
mid and low latitudes, the rhumbline is adequate for almost all calculations of 
course and distance. 





Outside this range, a more efficient track is the great circle, which is always 
the shortest route between two points on a sphere. However, in order to 
follow a great circle, a vehicle must be continuously changing its course. Since 
this is at best inconvenient, if not impossible, several rhumblines are often 
used to approximate a great circle. 





To plan a course using this technique, the navigator should first run the pro- 
gram Great Circle Plotting. For this program, the user must input the latitude 
and longitude of his starting point and his destination. Then, for any inter- 
mediate longitude A; that he specifies, the program will calculate the latitude | 
L; at which the great circle from source to destination will intersect the speci- | 
fied longitude. If several pairs of coordinates (Lj, ;) are calculated, then the 


next program, Rhumbline Navigation, may be used to find course and distance 
for the rhumblines linking these intermediate points along the great circle. 


The inputs to Rhumbline Navigation are the coordinates of two points on the 
globe; outputs are the rhumbline course and distance from the first point to 
the second point. The program may be used alone, to determine the rhumb- 
line from source to destination; or in conjunction with Great Circle Plotting, 
to compute several rhumblines to approximate a great circle. 


Intermediate points (Lj,Aj) 
_—_ tou by Great Circle Plotting 
Pa | 


<—_———Rhumblines——_—_—_> 






Great Circle__ 





(Ly 4 ) 


(L2,Aa) 








re 
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GREAT CIRCLE PLOTTING 


Equations: 


es tan L, sin (Aj —A,) — tan Ly sin (Aj - Az) z | 
oe sin Qa - A) 


where (L,,A,)= coordinates of starting point 
(L2, Az) = coordinates of destination 


(L;, Aj) = coordinates of intermediate point on great circle 


Note: 


The program does not compute along lines of longitude (A, = A). 
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sing = sin (y - Aa) 





tan, sin, 




















DEN =sin (\ -A, 


Display Lj in D.MS 
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a 


Key in program 





Input coordinates of starting 





point: 


Latitude ( ES for S) IDs L,, dec. deg. 
Longitude ( IQ for E) JD! Ai, dec, deg. 











Input coordinates of destination: 


Latitude ( for S) ,D. 








Longitude ( {iq for E) “iD: Az, dec. deg. 


Return to top of memory 








Input the intermediate longitude 





(EX for S) and compute cor- 





responding latitude 





For new intermediate longitude, 





go to step 5; for new source (or 





destination) go to step 2 (or 





step 3). 
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RHUMBLINE NAVIGATION 


Equations: 
m(Ay -X 
CH oe sa) ; 
180 [in tan (45 + Cy L,) - In tan (45 + = L,)] 
60 (Az -A,) cos L; cosC = 0 
D= 
(L2 -Li) 
cos C 


| 60 ——-———_ ; otherwise 
where (Ly, A;) = coordinates of initial point 

(Lz, Az) = coordinates of final point 
| C = rhumbline course 


D = rhumbline distance 
Notes: 
| 1. Nocourse should pass through either the south or north pole. 


2. The course may not go due east or due west across the 180° meridian 
(International Date Line). 


3. Errors in distance calculations may be encountered as C approaches 90° 
or 270°. 


4. Accuracy deteriorates for very short legs. 
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REGISTERS 





(dec. deg.) 











Let a= 1/2 (Ay -A,) {dec. deg.) 
Normalize a so that 
-180 <A; - Aq < 180; 
finds shortest route (dec. deg.) 
round earth 





























Let y = 10/90 

Let x =In tan, -In tan, 

C= tan"! y/x In tan 
(45+L3/2) 















































isplay course 
compute distance D 











If cos C #0, 
go to tine 45 

Cos C = 0; heading is 
due E or due W 




















O = 60 (Ay - Aa) cosh 
Halt and display Dist 
Heading is not due E or W 
Apply formula: 

D = 60(L, - £1 )/cos C 























ad ad | 


aad | 
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INPUT 
sree] wsrmcrone | oahtns ee ostons 


Key in program 





Input the initial latitude ( [CHS | 





for S) 














}3 | Input the initial longitude ( 
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Example: 


A ship sailing from San Francisco (L37° 49’N, \ 122° 25'W) to Tokyo 
(L35° 40’N, A 139° 45’E) will follow three rhumblines to approximate the 
great circle route. The navigator chooses the two intermediate points to be at 
 155°W and  175°E. Find the rhumbline courses the ship should follow, 
and the distance covered on each leg. 


Solution: 
First key in Great Circle Plotting. 


37.49 [5] PH] (0) 122.25 [3] Px) Bie) (2135.40 [5] PH) G2) 139.45 
[9] H) EGR [3] [4] LPram] 155 47.4606 


175 BS (es) ——_|@|_ —______—_— 47.3610 





Thus the two intermediate points are (LA7° 46'N, \ 155°W) and (L47° 36'N, 
X 175°E). 

















x | comers | 
Set flag to 0 for f(x) 


Recall x and branch to 
calculate f(x) 


Roll down to remove flag 














May Pause to see convergence 
Set flag to 1 for f(x + 8) 











| 2401 frcut_| 
EEX 

5 

$TO3 

+ 














ola} alalx] x) a) a 




















used to find f(x} and 
f(x +8), Flag in Ro is 
O for f(x), 1 for 


























f(x) /(x +) fs function value = 0? 


f(x}/Oc + 8) Yes, output solution 
No, check flag 











nN 
8 








~ 





4 





= 








SHR - 1) 





a 
8 














LUNE | 
00 
[ot | 
| 02 | 
03 
04 
0s 
06 
o7 
08 
08 
[10 | 
: 
[17 | 
[18 | 
| 19 | 
Hi 
| 26 | 
H 
= 
[38 | 
36 | 
37 
38 
38 
| 43 | 
= 
[a7 | 
| 48 | 
[49 | 
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OUTPUT 
DATA/UNITS 







INPUT 
step INSTRUCTIONS DATA/UNITS 
fa | Key in lines 1-16 of program 
[2 | Key in function f(x) 


Key in a branch to line 31 

























Press until display shows 





line 30 













Key in lines 31-49 of program 





Switch to RUN 













Store initial guess for solution 





Store tolerance 













Compute solution 









To change x; or € go to appro- 








priate step and store new value. 


Example: 


An equation often solved by gear designers is 
tanx-x-I=0 


where x is an angle in radians and I is the involute of x. Find the angle xo 
corresponding to an involute of 0.0324. 


Note: 


Since a gear designer might want to calculate x for several values of I, it will 
be simpler to store I in Rz for use by the function f(x). 
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Solution: 


Example User Instructions 


se] __mermuenons | outs ee 


Key in lines 1-16 of program 





Key in steps for f(x) = tan x - 




















Key in branch to 31 





Press ESid 8 times, until display 























Set tolerance € = 107° 





Compute solution xo 





Convert the angle to degrees 








Display last value of f(x) 


Last f(x) = 2.30 x 107° 
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NUMERICAL INTEGRATION, SIMPSON’S RULE 


Let Xo, X1, «5 X, be equally spaced points such that x; = x9 + ih for 
i= 0, 1, 2, ..., n at which corresponding values f(x), f(x;). ..., f(x,) of a 
function f(x) are known. This function need not be known explicitly but if 
it is, these values can be found previously by writing the function into 
memory and evaluating at the various points. n must be an even positive 
integer. 


Simpson's Rule is: 
p 


jf. f(x) dx = ; [f(xo ) + 4f(x1) + 2f(x2) +... + 4f(Kq_3) + 2F(Xp_2) 


xo 
+ 4f(Xn-1) + f(xn)]- 


Let the solution be indicated by I. 


ea 

a aa 
24 = RCL 1 

a 











ae [soa 


ze 01 | STOi 
R/S 


00 RCLO 

















24 7 RCL 1 


2301 | STO 1 


51 
23 _a o 1 
R/S 
24 = RCLO 























. 
Nn 


a a i 





| LINE | 
pa 
| 04 | 
05 
06 
07 
io" 
= 
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Key in program 





Store increment 


Enter first function value 
Enter last function value 
Enter values i = 1, 2,...,9-2 


Enter valuei=n-1 

















Example 


Tv 
Compute if sin? x dx using Simpson’s rule with h = 7/8. 


0 


The following data must be found first: 








Solution: 


iT 
iA sin? x dx = 1.5708 
0 


The exact solution is 7/2. 
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NUMERICAL SOLUTION TO DIFFERENTIAL EQUATIONS 


This program may be used to solve a wide variety of first order differential 
equations of the form 


y’ = f(x,y) 
with initial values Xo, Yo. 


The solution is a numerical solution which calculates y; for x; = Xo + ih, where 
h is an increment specified by the user andi=1,2,.... 


The program uses a modified Euler method (predictor - corrector): 
A _ h A 
Vier =¥i thf (x, yp Va [£j, vi) + (X41, Vir 1 


f(x, y) is keyed into memory starting at line 18. The user has 13 program 
steps to write f(x, y); registers Rs, Rg, and R; are also available. The user 
should assume that x and y will be in the X- and Y-registers, respectively. The 


routine should return with the value of f(x, y) in the X-register and should 
end with a GTO 31. 

















RCL 4 
1571 |gx-0 | 
1306 | GTO 06 
22 |RL 
2403 |RCL3 
+ 














oa 
= 





a 
oO 
rc 
o 
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Key in lines 1-17 of program 





Key in function f(x, y) 





Key in branch to line 31 





Press repeatedly until dis- 





play shows line 30 





Key in lines 31-49 of program 








Switch to RUN 


Store increment 











— initial conditions 


Display next x- value and cor- 








responding y- evalua 








Repeat step 9 as often as desired 


Example: 


Solve numerically the differential equation y’ = x \/y with initial conditions 
Xo = 1, Yo = 1. Use a step size of h=0.1. 





Solution: 


Key the function in as [ft] [x] 
ofa fa [is | a 
1.1077 | 1.2319 | 1.3745 | 1.5372 | 1.7221 
(by nen 





1.0 | 1.1078 | 1.2321 | 1.3748 | 1.5376 | 1.7227 
(ene’) 
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LINEAR INTERPOLATION 


If (x,, f(x,)) and (x2, f(x2)) are two points of a function f(x), then the 
function at Xo can be approximated by the following formula: 


(2 - Xo) f(x1) + (Ko - x1) f(x2) 


fia) = (x2 - 1) 


This is called the linear interpolation formula. Of course, x2 cannot equal x,. 


5 


iJ 
































2401 } RCL1 
61} x 


24 02 | RCL2 
2400 | RCLO 


~ 











ry 
“ 





.y 
we 





a) > 
&/e 





a 
a 





3 





=| > 
Orn 





& 





SUE ATT ET ATT 


D 
a 
| 03 | 
| 04 | 

05 
06 
o7 
| 08 | 
o9 
10 
re 
ae 
| 15 | 
16 
17 
18 
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INPUT 


Key in program 





Store first point 








Store second point 








Key in xo, find f(xo) 





Repeat step 5 for as many x- 





values as desired. 


Example: 
Given 
£(7.3) = 1.9879 
f(7.4) = 2.0015, 
find by linear interpolation f(7.37). 


Solution: 
(7.37) = 1.9974 


OUTPUT 
DATA/UNITS 





| : 


CHAPTER 6 STATISTICS 


CURVE FITTING—LINEAR REGRESSION 








When investigating the relationship between two variables in the real world, it 
is-a reasonable first step to make experimental observations of the system to 
| gather paired values of the variables, (x, y). The investigator might then ask 
the question: What mathematical formula best describes the relationship 
| between the variables x and y? His first guess will often be that the relation- 





ship is linear, i.e., that the form of the equation is y = a, x + a9, where a; and 
@p are constants. The purpose of this program is to find the constants a, and 
ao, which give the closest agreement between the experimental data and the 
equation y = a,;X + ag. The technique used is linear regression by the method 
of least squares. 





The user must input the paired values of data he has gathered, (xj, yi), 
i= 1, ..., n. When all data pairs have been input, the regression constants a; 
and ag may be calculated. A third value may also be found, the coefficient of 
determination, 17. The value of 1? will lie between 0 and 1 and will indicate 
how closely the equation fits the experimental data: the closer r? is to 1, the 
better the fit. 


Equations: 


y =a,X ta 


All summations below are performed fori=1, ..., n. 


Regression constants: 


| 





if- a .UU!”™””™S™CtOt:”t—<C—~—SSC( 
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where yooh 
—_ 2X 

x =— 

n 


Coefficient of determination: 


2 
E _ may 
he n 


T ee 
& 7 oe) & 7 oe) 
Note: 


The values for ag and a, are stored in Rg and Rj, respectively. After the 
calculation of ag, a,;, and r”, the estimated y-value, }, corresponding to any 
x-value may be calculated by y = a,x + a9. 





Programming Remarks: 


The intermediate value C = Yxy — (LxZy/n) is first calculated at line 14 but 
is also needed near the end of the program to find r?. Since all registers Ro 
through R, are in use, the only place to save this value is in the stack. Hence 
C is preserved in one or more of the stack registers from lines 14 through 36, 
when it is used. It is due to the presence of C in the stack that users are 
warned not to disturb the contents of the stack after calculation of ag and a; 
(see step 4 of User Instructions). 





MOA Wy 


Chapter 6 Statistics 89 


REGISTERS 


























Steps 1-7 for summation 





2405 | RCL5 
2407 | RCL7 



































C = Exy - (2x Zy/n) Es 























Halt to display ao 











Halt to display a, 


























E = Ey? - ((2y)?/n) 
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INSTRUCTIONS DATATUNTS 


Key in program 





Initialize 


Perform for i = 1,..., n: 





Input x-value and y-value 


Compute regression constants 
Compute coefficient of deter- 
mination 

To calculate a projected y-value, 
input the x-value 

Perform step 6 as many times as 
desired 

For a new case, go to step 2. 

* The contents of the stack 
should not be disturbed at these 
points. 















































Example: 


An eccentric professor of numerical analysis wakes up one morning and feels 
feverish. A search through his medicine cabinet reveals one oral thermometer 
which, unfortunately, is in degrees centigrade, a scale he is not familiar with. 
As he stares disconsolately out his window, he spies the outdoor thermometer 
affixed to the windowframe. This thermometer, however, will not fit com- 
fortably into his mouth. Still, with some ingenuity.... 


The professor suspects that the relationship is F = a, C + ag. If he can geta 
few data pairs for F and C, he can run a linear regression program to find a; 
and ao, then convert any reading in °C to °F through the equation. So tossing 
both thermometers into a sink of lukewarm water, he reads the following 
pairs of temperatures as the water cools: 


100 


If the relationship is indeed F = a,C + a9, what are the values for a, and ag? 
What is the coefficient of determination? 

















[102 | 


102 
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Solution: 


(f] Ceram] [1] 40.5 [4] 104.5 [Rs] ———————» _ 1.00 
38.6 [4] 102 [Rs] ———_—_——______"___» 2.00 
37.9 [4] 100 [ws] —————_—_ 3.00 
36.2 [4] 97.5 [Rs] ————--______——>_ 4.00 
35.1 [4] 95.5 [Rs] ———————————__ 5.00 
34.6 4) 94 [Rs] -—_——__—_-_——————>_ 6.00 
Gq (9) (8) (Ss) —_- > _ 33.53 
[RIS] 
es | 0.99 


Thus, by the data above, F = 1.76 C + 33.53, with r? = 0.99. (The real equa- 
tion, of course, is F = 1.8C + 32.) 


Suppose the professor puts the centigrade thermometer in his mouth and 
finds he has a temperature of 37°C. Should he be worried? 


37 GSS 2) (x) GSS (o] |] ————————————_>_ 98.65°F 


It looks like he is safe. 
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EXPONENTIAL CURVE FIT 


This program computes the least squares fit of n pairs of data points {(x;, yi) 
i=1,2,..., n}, where y; > 0, for an exponential function of the form 


y=ae* (a>0). 
The equation is linearized into 
In y =Ina+ bx. 


The following statistics are computed: 


1. Coefficients a, b 


1 P 
2x; In yj - a (2x;)(Z In yj) 





Dx;? - - (xj) 








2. Coefficient of determination 


2 


bx In yj - a Lx; Z In 7 


bs _ Gx) | E dnyy? ny | 
n n 


3. Estimated value f for a given x 





J=ae* 


Note: 


nis a positive integer andn #1. 
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DISPLAY KEY DISPLAY | reasons 
rae [oe] = | 
Roa 
: 


| 00 | 

jot | 

poz | 31 | 
| 03 | 1502 | 

| 04 | 2351 02 | 

| 05 | 

| 06 | 

| 97 | 





R,Ziny 
Rs IxiIny 


ey 


= 
a 
oO 
~“ 


24 07 





| 2406 | RCLE | 


71 
41 
18 1502 | gx? 
19 2403 | RCL3 
2 14+ 
71/4 + 


pa {atl 


| 24 | 2so7{ rci7_| 

















OUTPUT 
DATA/UNITS 


Key in program 





Initialize 


Perform for i = 1,...,n: 





Input x-value and y-value 





Compute constants 








Compute coefficient of deter- 





mination 





To calculate ¥, input x 








Perform step 6 as many times 





as desired 





For new case, go to step 2. 








* The stack must be maintained 


lamin 
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Example: 


fs | 72 | 131 | 9s | 258 | 3.14 | 
fyi | 216 | 161 | 116 | 85 | 05 | 


Solution: 
a=3.45,b=-0.58 

y = 3.45 e70-58x 

1? =0.98 

For x= 1.5, = 1.44 
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LOGARITHMIC CURVE FIT 


This program fits a logarithmic curve 
y=at+blInx 
to a set of data points 
{i yi, i= 12553 n} 


where x; > 0. 


Program computes: 


1. Regression coefficients 


Syuliie a hing So: 
n 
Wie atom Caen eer he 
>> (In x;)* =e (z In x)" 


a= — (Zyj-bEInx) 


2. Coefficient of determination 


2 
[Py inxi-22In x Dy 


r= 


3. Estimated value ¥ for given x 


y=atblinx 


Note: 
nis a positive integer andn # 1. 


95 














96  Chapter6 Statistics 


DISPLAY KEY 


| CODE | ENTRY 
MMA 


31 



























REGISTERS 
Roa 
R,b 








Rs Ly Inx 


























Key in program 








INSTRUCTIONS 









Initialize 


Perform for i= 1,...,.: 





a 





Input x-value and y-value 








Compute constants 








Compute coefficient of deter- 





mination 





a. 
To calculate y, input x 








Perform step 6 as many times as 





desired 





For new case, go to step 2. 








* The stack must be maintained 





at these points 












R, Ly? 
4 
5 
6 


Re Zinx 


OUTPUT 
DATA/UNITS 
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Example: 


pm| 3 | 4 | 6 | 10 | i | 
| vijis | 93 | 234 | 45.8 | 60. | 


Solution: 

a =~47.02,b = 41.39 

y = -47.02 + 41.39 Inx 
17 =0.98 

For x = 8, ¥ = 39.06 
For x = 14.5, 9 = 63.67 
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POWER CURVE FIT 


This program fits a power curve 


y=ax? (a>0) 


to a set of data points 


{(i, yD, i= 1, 2, sath 


where x; > 0, y; > 0. 


By writing this equation as 


Iny=bInx+tlIna 


the problem can be solved as a linear regression problem. 


Output statistics are: 


1: 


Regression coefficients 





‘ ] : 
Suiting OO 
b= n 
2 
> (In xj)? - ae 


2 In yj 2 In x; 
a= exp =o ee 


Coefficient of determination 


2 
(Z In x;) (2 In S 


n 


s (In x;) (In yj) - 
2 = 





r= 


n 


E as. 2a E fii ge ee 
1 n 1 


Estimated value y for given x 


<> 
! 
rs) 
*< 


Note: 


n is a positive integer and n# 1. 








~~ 6 tee $00 08 Foot SR as 5S eee ee et ee a ee ater cai i, ————— 


PLAY 


24 07 A aa aed 





R, = (In y)? 





Rs = (In x) (Iny) 


R,DInx 





~ 
= 





= 
a 
Oo 
~ 





4 


n wl+|x 

fe) i) 
° ae cr 
° 7) > 


Dip 
~ 
un 





2404 | RCL4 
61 |x 


| 2403 | RCL3_| 


3\" |Z\8 8 
Pad 

124 Dar, 
nN 


~ 
= 


~ 
= 


| 2401 | 
at | 
i Be | 
= 8h | 
| 2402 | 
| 2404 |RCL4 
| 1502 | 
| 2403 | 
ie 
Pt | 
ee 
| 1300 | 
fee = — 1 
Sa 


INPUT OUTPUT 
STEP INSTRUCTIONS DATA/UNITS DATA/UNITS 








Input x-value and y-value 





Compute constants 








Compute coefficient of deter- 





mination 








Input x-value and compute ] 








Perform step 6 as many times as 





desired 





For new case, go to step 2. 








“ The stack must be maintained 








at these points. 
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Example: 


fx; ] 10 [ 12] 15] 17[ 20] 22 | 25 30 | 32 | 35 
1.25| 1.41| 1.73| 2.00 |2.53] 2.98} 3.85 | 4.59 | 6.02 


Solution: 

a= .03,b= 1.46 

y= 03x!46 

r? =0.94 

For x = 18, ¥ = 1.76 
x=23,9 =2.52 
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COVARIANCE AND CORRELATION COEFFICIENT 


For a set of given data points {(i, yj), 1= 1,2,..., n}, the covariance and the 
correlation coefficent are defined as: 


covariance Sx y =! (o - ~ Lx,Lyj 


r_| | 
OF Sxy = fe os >) 


‘ Sx 
correlation coefficient r =—~ 
Sx Sy 


where s, and sy are standard deviations 


Dx;? - (2x;)?/n Lyi? - (Zyj)?/n 
eT n=1 Sy ~ n-| 
Ron-t 


R, ry? 





So 
wo 
: 





+ 


ever ae 
+ Z 
= Z 





Z 
g 
Z 


Q 
fs) 
8 





2404 | RCL4 
CL7 
Bt | 
| 2403 | 


B:) 


Alri an Pl a 
a oO oO 

Cc r 
= wo wo 





r 
ee c 
z 


TOO 
71 | + 
a0 | 


eae RCL 3 





x|/vin 
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2 Sh eee ee 
Key in program —— Za a as ——— 
Initialize Se eee 


Perform this step for i = 1, 2,...,n 








Compute covariance Sxy 





Compute sxy’ 





| Compute correlation coefficient 





For new case, go to step 2. 


Example: 

[xi | 26 | 30 | 44 | 50 | 62 | 68 | 74 | 
[vil 92 | 8s | 78 | 81 | 54 | 51 | 40 | 
Solution: 

Sxy = -354.14 


as = -303.55 
r= -0.96 
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MOMENTS AND SKEWNESS 


This program computes the following statistics for a set of given data {x1 »X2;5 
ies Xn}: 


18* moment 


. 


1 = 
2"4 moment m,=— 2x;? -x? 
n 


1 3 = 
34 moment — m3 = a Dx;? - — X Lx,” + 2x? 
n 


moment coefficient of skewness 


Di DI 


(77) 
Uv 
e 
> 
=< 


= 
z 
m 





2407 | RCL7 
2403 | RCL3 


cas 
jax? 
ee 
23.01 


2403 | RCL3 
ri cs 
2407 | RCL7 


Cae 
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a in Program 
Initialize = 





Perform for i = 1, 2,.. 





Input x-value 


Delete erroneous data 








Compute the mean 





Compute the second and third 





moments 








Compute the moment coefficient 





of skewness 





For new case, go to step 2. 


Fifi} 2{3]4[s]o6]7 Ja lo | 
[xi [2.1] 3.5] 42] 65] 41] 3.6] 5.33.7] 49 | 


Solution: 
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NORMAL DISTRIBUTION 


The density function for a standard normal variable is 


‘je es, 


Vin 


The upper tail area is 


> 


ZZ 


For x > 0, polynomial approximation is used to compute Q(x): 
Q(x) = f(x) (b; t+ bz t? +b, t? + by t* + bs t*) + e(x) 


where (€(x)| << 7.5 x 1078 


po. r=6.39164i9 
[ee tx 


b, = 31938153, b2 = -.356563782 
b3 = 1.781477937, bg = -1.821255978 
bs = 1.330274429 

Note: 


The program only works for x > 0. Equations f(-x) = f(x), Q(-x) = 1-Q(x), 
where x 2 0, can be used to find f and Q for negative numbers. 


Reference: 


Abramowitz and Stegun, Handbook of Mathematical Functions, National 
Bureau of Standards, 1968. 
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DISPLAY DISPLAY KEY 








STO7 
R/S 


RCL O 
RCL6 
x 
1 
+ 


g 1/x 

a 
St 
31 | f 


INPUT OUTPUT 
STEP INSTRUCTIONS DATA/UNITS KEYS DATA/UNITS 


Key in program 
































00 | 
Por | 
02 | 
08 | 
04 | 
Fo5 | 
06 | 
For | 
08 | 
09 | 
10 | 
Fatg| 
rae | 
Ae 
18 
19 
2a | 





Initialize 





Store constants 











Input x and compute f(x) 





Compute Q{x) 





For a new case, go to 4. 
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Examples: 
1. x=1.18 


2. x=2.28 
Solutions: 


1. = f(x) = 0.20 
Q(x) = 0.12 


f(x) = 0.03 
Q(x) = 0.01 
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INVERSE NORMAL INTEGRAL 


This program determines the value of x such that 


t? 
ae § 
a- f ait 
x 20 


where Q is given andO0<Q<05. 


The following rational approximation is used: 


Co tC, ttc, t? 
x =t -——________——_ + &(Q) 
l+d,tt+d, t? +4, ¢° 

where {e(Q)| <4.5 x 1074 

1 
n — 

Q? 
Cg = 2.515517 = d, = 1.432788 
c, = 0.802853 = d, = 0.189269 
¢2 = 0.010328  d; = 0.001308 


t= ! 


Reference: 


Abramowitz and Stegun, Handbook of Mathematical Functions, National 
Bureau of Standards, 1968. 
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DISPLAY KEY 
ENTRY 























2400 | RCLO 
RCL7 


2407 












= = 


: 


M\wjiwrmlo |x 


x >|+|> 
ro 4 
c io} 
a a 






Ry 1+ditt+dt? tds? 

















= 3 


BE 
| 61 | 
| ot 
aor | 7 
CLX 
a 02 | RCL2 


eed CLI 




















Bue see el eee ee 


OUTPUT 
DATA/UNITS DATA/UNITS 


Key in program 





Initialize 





Store constants 


























ls | For a new case, go to 4. 
Examples: 

1. Q=0.12 

2. Q=0.05 


Solutions: 
1 = x=1.18 
2. x=1.65 
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FACTORIAL 


This program will compute factorials for positive integers between 2 and 69. 


n! =n(n-1)(n-2)...(2) (1) 
Notes: 


1. For large values of n, the program will take some time to arrive at a 
result, up to a maximum of about 20 seconds for n = 69. 


The program does not check input values and will return incorrect 
answers for values of n < 2 or n > 69 or n non-integer. 


DI 























RCLO 
GTO 00 











| LINE} CODE | 
|_oo f 2] 
ae | 26] 
Lor be 
| 03 | 28 
04 29 
os | 30 
06 
07 
08 
o9 
| 10 | (ap 
36 
ea =} — 
13} | 38] 
pe es} 
ca ao 
| 18 | 
| 19 | | aaf 
| 20 | 45] 
|_| | 46{ 
Ea a7 {| 
| 23 | | a8] 
| 24 | El 


24 00 
—— 
ae 
=a 
ce 
se 
hs ed 
we 
el 
ee al 





a ee ee ee eee ee 


Initialize 








Key inn (2<n<69) 





For a new n, go to step 3. 


Examples: 
1. 5!=120.00 


2. 10! = 3628800.00 


Chapter 6 Statistics 


Ill 
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PERMUTATION 


A permutation is an ordered subset of a set of distinct objects. The number of 
possible permutations, each containing n objects, that can be formed from a 
collection of m distinct objects is given by 


=m(m - 1)...(m-n+t 1) 


7 ~ pl 





where m, n are integers andO<n<m. 


Notes: 


1. .,P, can also be denoted by Pf’ , P(m,n) or (m),. 


2. mPo =1,mP; =m, mPm =m! 


DISPLAY KEY 
ENTRY 


1315 | GTO 15 


R 


GTO 00 


GTO 00 
1 


REGISTERS 


[Rom 
PT 
i 
Rs 
Ra 
Rs 
Re 
R; 





= 








nN 
nN 
<_ 





o 
= 








41 


1871 |gx=-0 _| 
13 37 


rc 
nw i 
2 


13 39 | GTO 39 
0 
































f LASTx 





RCLO 
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INPUT OUTPUT 
mermucrons | oats | MS st | 


Key in program == ae L JL | 





Store m,n 








Compute permutations 





For new case, go to step 2. 


Examples: 
i. 43 P3 = 74046.00 


2. 73 P4 = 26122320.00 





—_—_. ALP 
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COMBINATION 


A combination is a selection of one or more of a set of distinct objects 
without regard to order. The number of possible combinations, each 
containing n objects, that can be formed from a collection of m distinct 
objects is given by 


(<= m! _ m(m-1)...(m-n+ 1) 
mn (m =n)! n! } 2°" 4 


where m, n are integers andO0<n<m. 





This program computes ,,C,, using the following algorithm: 


1. Ifn<m-n 


Cc arnt! ment? . mm | 
mCn 1 ge es | 
2. Ifn>m-—n, program computes mCm_n- 


| Notes: | 


1. mC,, which is also called the binomial coefficient, can be denoted by | 


CM , C(m,n), or (f). 





2.0 mCn =mCm-n 











= oe a oe ee 
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[cove | &NTr 


1473 | f LASTx 
1441 | fx<y 


| | 
ol 
| 2301 | s 


51 












































+ 


ae 

2307 [S70 
xzy 
fx2y 

| 13.22 | 

| 2402 | 








TO 22 
RCL 2 
GTO 00 
| 22 | 
| 2400 | 
| sit | + 





| 00 | 
| on | 
| 02 | 
| 93 | 
as 
| 07 | 
08 
09 
10 
11 
| 12 | 
| 13 | 
ie 
| 19 | 
| 20 | 
| 21 | 
[22 | 
| 38 | 
| 24 | 


INPUT OUTPUT 


For new case, go to step 2. == 


Examples: 
1. 43C4 = 1088430.00 


2% 4303 = 12341.00 


xz 
®D 
< 
5 
3 
8 
3 
3 
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CHI-SQUARE EVALUATION 


This program calculates the value of the x? statistic for the goodness of fit 
test by the equation 


(O; - E\)? 
Ej 


where 0, = observed frequency 


E; = expected frequency. 


The x? statistic measures the closeness of the agreement between the 
observed frequencies and expected frequencies. 


Notes: 


1. In order to apply this test to a set of given data, it may be necessary to 
combine some classes to make sure that each expected frequency is not 
too small (say, not less than 5). 


If the expected frequencies E; are all equal to some value E, then E 
should be computed beforehand as 





and then input at each step as the expected frequency E;. 












DISPLAY 


=] 





Zam: 
Zue 
Z 
Z 
e 
Z 
Z 


2300 | STOO 
2301 | STO1 








‘oO 
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ilala 
ala 
nN 


2402 | RCL2 


o 














= 


O | 235101 | STO+1 








TO 04 


TO~ 1 
CLO 


~ o o|o|o 5 
“ o NO) a mi 


LINE: | 
25 | 
| 26 | 
| 28 | 
| 29 | 
| 30 | 
| 31 | 
= 
| 37 | 
| 38 
| 39 | 
| 99 | 
aa 
| 42 | 
| 43 | 
| 44 | 
| 45 | 
| 47 | 
| 48 | 
| 49 | 





= = Bb 


Din z n o]+ 
ro) 4(S)S 3 
r ° ° e 
nN NR o o 





STEP INSTRUCTIONS 


INPUT 
DATA/UNITS 


Key in program 





Initialize 


Perform for i = 1,..., 1n: 





Input observed and expected 





frequencies 





Delete erroneous data 








Display x? 





For new case, go to step 2. 
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Solution: 


x? = 4.84 
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PAIRED t STATISTIC 


Given a set of paired observations from two normal populations with means 
M1, M2 (unknown) 


The test statistic 


which has n - 1 degrees of freedom (df), can be used to test the null 
hypothesis 


Ho: Ha = be. 
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REGISTERS 
R 2 
GTO 00 


c 

Zz 

m 
3 
< 


y 


Z 
Z 
A 
Z 
Ae. 
ZEEE 
Z 














z+ 
s 
Bea ae 
| taaitee | 
ptf ey | 
1 
GT 








7 
74 {| R/S 





1 
2403 | RCL 3 
01 

1300 | GTO 00 
a eee 
aa Se 























| 00 | 
roe 
| 04 | 
| 05 | 
| 06 | 
| 08 | 
| 09 | 
| 10 | 
oa 
Le] 
| 15 | 
| 16 | 
= 


[a|2|3]a|a|a|e|s||e|s 


INPUT OUTPUT 
Key in program 
Initialize 


Perform for i = 1,..., n: 
Input one pair of observations 
4 


Delete erroneous data 

















Compute t and df 








For new case, go to step 2. 
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Example: 


js [4 | 17s | i7_| 17s | 154 


Solution: 
t= -7.16 
df = 4.00 





nee nn en 
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t STATISTIC FOR TWO MEANS 


Suppose {x1, X95 ics Xn,} and {y15 Ve ca Yno} are independent random 
samples from two normal populations having means f,, “2 (unknown) and 
the same unknown variance o?. 


We want to test the null hypothesis 
Ho: #1 - 42 =D 
where D is a given number. 


Define 





X= Xj 
ab i=1 
we 
y=-— Dy 
M2 j=1 
fu x-y-D 
bs Ex? -n, X? + Zy;? -n, y? 
Nn, ng ny +n -2 





We can use this t statistic, which has the t distribution with nj +n, ~- 2 
degrees of freedom, to test the null hypothesis Hg. 
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DI 


a 


PLAY DISPLAY 
ENTRY LINE 
2403 | RCL3 2402 | RCL 2 
2406 | RCL6 
142 


STO 3 
STO6 
23 07 


/S 


Z 
g 
Z 




















+ 


gZ 
Z 
: Z 
£ = Z 
ti] + |] jpDIe |D RIBI+ 1 AI7 12 
S| |slPisle| 
w ° 


> 
“ 
Ny 
~ ~ o|s 
= —_ o|= 


LINE 
roe 
05 | 
| 96 | 
or | 

08 
[10 | 
Es 
12 
13 | 
ace 
15 | 
16 | 
| 47 | 
18 | 
| 19 | 
| 20 | 
| 2 | 
| 22 | 
oe 
[24 | 


71 


ec 
Key in program ——Ee 
Initialize [ij | = 

ae as ee | See 


re eer’ 


Perform for i = 1,...,: 





Input x-value Xj 


— irons J 


5 | Perform for i= 1,..., nz: ee he 
Input y-value 2 es a | 
a i en ea a 


Input D and compute t 








To find the means of x- and y- 





values 














For a new case, go to step 2. 





126 Chapter6 Statistics 


Example: 
x: 79, 84, 108, 114, 120, 103, 122, 120 
91, 103, 90, 113, 108, 87, 100, 80, 99, 54 
n, =8 
n, = 10 
D=0 (ie., Ho: Ma = M2) 


Solution: 
t= 1.73 


x = 106.25 
y =92.50 
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ONE SAMPLE TEST STATISTICS FOR THE MEAN 


For a normal population (x;, X2 ..., X,;) With a known variance o”, a test of 
the null hypothesis 


Ho: mean p= Uo 
is based on the z statistic (which has a standard normal distribution) 
vn (X - Uo) 
z= ——_—__—_ 
oO 
If the variance o? is unknown, then 
= va (X = Ho) 
s 


is used instead. This t statistic has the t distribution with n - 1 degrees of 
freedom. X and s are the sample mean and standard deviation. 


ISPLAY DISPLAY 
REGISTERS 


DI 


NZ 
PC 
Z 


x 
N 
< 


[Rs Used 


Z 
reales EAR 
~ x Z 


~s 
slala)*|lalalelg|*|zlz 
< mi?) x r 
Oo w 
= 
- a > 
z 


+ 


LINE 
| 00 | 
| on | 
| 02 | 
| 03 | 
| 04 | 
| 05 | 
| 06 | 
ae 
| 08 | 
| 09 | 
[10 | 
1 | 
| 12 | 
|_| 
aa 
| 15 | 
| 16 | 
17 | 
| 18 | 
| 19 | 
| 20 | 
| 24 
| 22 | 
mk 
|_24 | 


= 
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INSTRUCTIONS Dare niTe Pes ae 


Key in program 

Initialize = 
Perform for i = 1,. 

Input value 


Input Wo 











Compute t 





or 


Input o and compute z 





For new case, go to step 2. 


Example: 
Suppose Uo = 2, for the following set of data 


{2.73, 0.45, 2.52, 1.19, 3.51, 2.75, 1.79, 1.83, 1, 0.87, 1.9, 1.62, 1.74, 1.92, 
1.24, 2.68 } 


Solution: 
test statistic t = ~.69 


or z= -.57 ifo=1. 
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CHAPTER 7 SURVEYING 


FIELD ANGLE TRAVERSE 


A traverse is a series of line segments joined with specific lengths and angular 
relations to each other. With many applications in surveying, the field angle 
traverse may be used in establishing boundary lines, road layout, and in 
numerous construction situations. The transit and “chain” (commonly a 
length of steel tape) are often used to establish the angles and distances 
involved in a field angle traverse. 


Starting at a known point from a given reference direction, the transit man 
establishes the direction of a new line by measuring the angle or deflection 
turned to align the scope of the transit to the new line. With a measured dis- 
tance to the end point on the new line and its direction, coordinates of the 
end point relative to the origin may be established. The transit is then moved 
to the new “origin”, the reference direction is the line just determined, and 
the process continues. 


To run this program, the user must input the northing and easting of his start- 
ing point, the reference azimuth, and then the direction and distance from 
each point in the traverse to the next point. The direction may be input either 
as a deflection right or left, or as an angle right or left. The distance may be 
input either as horizontal distance, or as slope distance with zenith angle. 


Equations: 

H Dist = S Dist sin (Znth ang) 

Ni+1 = Nj + H Dist cos Az 

Ej+1 = Ej + H Dist sin Az 

Area = 4[(N2 +.N,) (E2 - E,) + (N3 + No) (E3 - E2) + 
... +(N, +N.) (Ei -E,)] 








130 Chapter7 Surveying 


where: N, E = Northing, easting of a point 
Subscript i refers to current point 
Subscript n refers to next to last point 
Numeric subscript refers to point number 
Az = Azimuth of a course 
H Dist = Horizontal distance 
S Dist = Slope distance 
Znth ang = Zenith angle 


Notes: 


1. The calculation for area of a closed traverse may be inaccurate for cases 
in which the coordinates of the figure are quite large, such as in state 
plane coordinate systems. In such cases, the user may run the Area By 
Double Meridian Distance program to calculate an accurate value for 
area once the bearings and distances have been established by this 
program. 


All angular inputs and outputs are in the form degrees, minutes, and 
seconds (DMS). 
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180 + Az 





180 + Az 


Current 





Ni 


E 





Ni 


Initialize “previous N’” 





180 + Az 


Clear R3, Ra, for 


a ZH Dist 





180 + Az 
180 + Az 


accumulation 





180 + Az 





Convert to decimal degrees 


























Deflection comes in here 














+ 


* 
¥ 


x 








1409 


= 





a 
°o 
a 





g 

o 

+ 

Nn 
| 


‘il 





35104 |STO+4 





2401 


~ 
3 
— 
i 
ing 
> 





18 
20 
| 21 | 
Fae] 
|_25 | 
| 26 | 
|_27 | 
| 28 | 
| 29 | 
|_30 | 
| 31 | 
| 32 | 
| 33 | 
| 34 | 
| 35 | 
| 36 | 
| 37 | 
| 38 _| 
[39 | 
| 40 
| 41 | 
Par] 
[45 | 
| 46 | 
| 47 | 


qalais 
ole] %Jer 
apes} | as 
Nn 


AA = (Nj + Ning AE 


Display Northing 


Display Easting _ 
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INPUT 
sre] _memerons [fim | ame | nt 


Key i in program 





Input the starting point coord- 





inates 








Input the reference azimuth 





If angle right 





If angle left 





If deflection right 





If deflection left 





If horizontal distance 








. If slope distance, input Zenith 





Angle and Slope Distance 














Repeat steps 4-5 for successive 




















Example: 


The diagram below shows measurements taken for a closed traverse. Find the 
coordinates of points 2, 3, and 4, the total horizontal distance traversed, and 
the area of the figure. 






OL 100°24'55" 


ee a Pee et 
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Solution: 


150 ESQ) (1) 400 ESS) [2] 311.3955 [+] Ceram ] [Rs] ———+ 0.00 
113.3455 177.966 [R/S] —————_______» 224.515 (N;) 





561.615 (E) 
100.2455 GS (1) [9] 86.0139 [4] 161.880 


[2] [6] [p/s] 





356.526 (N3) 
—_—___________—__» 468.601 (E;) 
87.3559 203.690 > 232.335 (Na) 
————_____—___—_—\— 307.151 (E,) 
100.4559 ES) Ke) [7] 124.0 [Rs] ———————- 149.903 (N,) 
399.784 (E;) 
[3] ————_—————__» 667.144 (2H 
Dist) 
[4] ——__—_———————> 26590.68 
(Area) 









356.526 








232.335 Area = 26,591 Sq. ft. @ 
aor 1S 224.515 
® 561.615 


i : = 149.903 
Calculated ending coordinates -N = 399.784 
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AREA BY DOUBLE MERIDIAN DISTANCE 











This program computes the area of a straight-sided closed figure from the 
bearings and lengths of its sides. It is generally more accurate than methods 
which calculate area from the coordinates of the figure. 


Area = > Lu DMD, x Latitude; 


DMD; = DMD;_, + Departure;_,; + Departure; 


where 


Departure; = Dist; sin Az; Latitude; = Dist; cos Az; 


Note: 


Angles are input as bearing and quadrant code. The quadrant code is 1 for NE, 
2 for SE, 3 for SW, and 4 for NW. 


































































LINE | CODE | |UNE | CODE _| sient 
ee eee 
01 26 f>R R, DMDi-4 
aw xz R 2 Departure}. 
ex | 28 | 2402 |rci2 Rs Area 
p29 | at [xy 
| 05 | | 30 | 2302 |sto2 
06 1401 | fINT + 
07 1461 | f x#y = RCL 1 
13 165] GTO 16. es + 
| 34 | 2301 |stO1 
ee eee 
(30 .| oz 2 
Ea oe 
| 38 [235103 |sto+3 | 
ef see RCL 3 
| 13.00 | GTO 00 
es as 
pa {| 
7 Lan |e a 
zoo{stoo || 4{[ | | 
| 22 | 1400 | f>H.ms | | 47 haz 
eae = 
2400 | RCLO 
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Key in program 





Initialize 





Input bearing 





Input quadrant code 





Input distance 





Repeat steps 3,4,5 for succes- 





sive courses. Area is displayed 





after last distance has been input. 


Example: 


Compute the area of the figure below. 


Solution: 
Area = 20937.44 sq. ft. 


YJ Wiz0osizam ERR: 
CloizitSevyiorn Ss: 
Line O2-9 inseey ENTER 


Line Oe OG GWANGE To 4GTo ls” 
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INVERSE FROM COORDINATES 


This program uses coordinates to calculate distance and bearing between 
points of a traverse. The area in square feet and a summation of distance 
inversed are also computed. 


Epa By 
H Dist = V(N; - Ni_1)* + (Ej - Ej-1)? ising! — 1 
Ni - Ni-1 


Area = 4[(N2 + N,)(E2. - E:) + (N3 + Nz) (E3 - E2) + 
... (Nn + Ni) Ei -E,)] 


where N,E= Northing, easting of a point 
Subscript i referes to current point 
Subscript n refers to next to last point 
Numeric subscript refers to point number 
H Dist = Horizontal distance 
Az = Azimuth of a course 


a 
3 
e 
> 
< 


REGISTERS 


TO 2 





R2 Current E 
Rg ZH Dist 
4 


ENTRY 
= 
eee 
ae a 

74 | R/S + 
f SIN 
gSiN"' | 
g x<0 
CHS 











| 2a | x2v 


RCL 1 


41 |- 
235101 |sTO+1 


1809 |g>P Cl 1 
235103 |STO+3 TOO 


21 











[LINE | 
|_00 | 
jor | 
| 02 | 
| 03 | 
| 04 | 
95 | 
[96 | 
ca 
| 08 | 
| 09 | 
18] 
La 
12 | 
ae 
18 














R 


1306 |GTO 06 





nN 
NO 











i] 
nN 
we 
E 
= 
o|~m 
Bis 
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(INPUT OUTPUT 
marmcrons oats |S ont 


Key in program 





Input starting coordinates 








Input next coordinates and 





display distance 








Compute bearing and quadrant 





code 








Repeat steps 3-4 for successive 





courses 





Display total distance inversed 





Display area of closed figure 





(ignore the sign) 





Example: 


350.000 
325.000 





N 7°07'30" E 
201,556 








Begin N oes © 


Area = 20937.5 Sq. ft. 
Total distance inversed = 641.033 
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CHAPTER 8 
TRIGONOMETRY AND ANALYTICAL GEOMETRY 


COORDINATE TRANSLATION AND ROTATION 


There are occasions, for instance in cartography or metalworking, when it is 
necessary or advantageous to shift one’s frame of reference. In mathematical 
terms, the occasion calls for a translation and/or rotation of the coordinate 
system. The origin is translated from (0, 0) to a new point, (xo, Yo), and the 
x and y axes are then rotated through an angle a to give new axes, x’ and y’. 
Suppose that a point P has coordinates (x, y) with respect to the old system 
of x and y axes. The problem then is to find the coordinates (x’, y’) of P 
with respect to the new system whose axes are x’ and y’. The diagram below 
illustrates this situation. 











(X09, Yo) 








Equations: 


x’ =(X-X9) cosa+(y - yo) sina 


y' =x - Xo) sina + (y - yo) cosa 


| 
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Notes: 


1. The program may be used to solve a problem of translation only, or of 
rotation only, or of combined translation and rotation. If the problem 
involves translation alone, a value of a =0 must be input. For rotation 
alone, the values x9 = yo = 0 must be input. 


The program assumes the following sign convention: a should be input 
as a positive number if the rotation is counterclockwise, and negative 
if clockwise. 


Programming Remarks: 


This program demonstrates a particularly powerful application of the polar- 
to-rectangular conversion ( [4] [+R] ) when combined with the capabilities of 
the four-register stack. The subterms (x — xo) cos a, (x - xX) sin a, (y - yo) 
cos a, and (y ~ yo) sin @ are all generated through [+f] [>R] and stored in the 
stack until needed. A more straightforward program using [#] and [f] 
would have required 30 program steps (as compared to 19) and one 
more storage register. 
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DISPLAY 


Est2 tse 














REGISTERS 





























2403 


Ax sing 





2401 





Ax cos & 














Ee cre = 
4x cos @ Ay cos @ 








Ax sina 











= ~Ax sina + Ay cosa 
































ro 
00 
ot | 
por 
| 98 | 
06 
o7 
08 
09 
10 
Fs 
-e 
|_21 | 
|_22 | 
|_ 29 | 
|_24 | 
|_28 | 
|_26 | 
|_27 | 
| 28 | 
|_29 | 
[30 | 
| 3 | 
[32 | 
| 33 | 
| 34 | 
[35 _| 
Pe 
| 39 | 
| 40 | 
[44 | 
[42 | 
| 43 | 
4] 
| 48 | 
| 46 | 
| 47 | 
| 4 | 
|_40 | 








Ax = x = x9 
Ay =y-Yo 
, 
, 


XIX |b |b I> 
i<ji< |< 
alalo 
8/812 
e|e ir 

CS 
< 
2 
g 
R 
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INPUT fe) 


Key in program 





Store origin of new coordinate 





system 











Store angle of rotation 


Convert coordinates from old to 





new system 














Perform step 4 for as many 








points as desired 








For a new case, go to step 2. 


Example: 


A backpacker’s route will take him cross-country away from the marked trails 
of an area. He knows that he will have to check his compass frequently against 
his map over this terrain, and regrets that the map is in such an inconvenient 
format for his purposes. In the first place, the grid lines on his map represent 
distances in feet from an origin about 25 miles away, which are such large 
numbers that they are hard to calculate with. Secondly, the map’s grid is 
based on true north while his compass readings are relative to magnetic north, 
a variation of 17°. 


Before he leaves home, the packer decides to draw a rough version of the map 
for his own convenience, locating his origin at the grid point (54000, 
118 000) and rotating his axes by 17° in a clockwise direction. As a first step, 
he wants to find the new coordinates of the bridge and the peak of the hill, 
whose coordinates in the old system are (55 750, 119 300) and (57 450, 
120 500) respectively. 
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RSE 



















iran 
a 





SS 
Tp 
5, 













119 000 KOSS 
_RSISNCSES | 
54 000 55 000 56 000 ee 57 000 
True Mag 
7 ' 
Solution: 


54000 Eig [0] 118000 Eig) [7] 17 Lets) Eiey (2) [+] [Prom | 
55750 [4] 119300 a 1293.45 
1754.85 


The new coordinates of the bridge are (1293, 1755). 








57450 [4] 120500 [R&S] ——————______» 2568.32 
__— ee YT 





The new coordinates of the peak are (2568, 3399). 


| 
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TRIANGLE SOLUTION B, b, c 


Given two sides and a non-included angle, this program solves the triangle for 
the remaining parameters by the following formulas: 


in B 
i. c= oe ) 
sin! (| —— 
2. A=2sin"' 1 -(B +C)=7 radians - (B + C) = 180° -(B+C) 


= 200 grads - (B+ C) 
- bsin A 
sin B 


If B is acute (<< 90°) and b <c, a second set of solutions exists and is 
calculated by the following formulas: 


4. C’'=2sin7' 1-C 
5. A’=2sin7! 1-(B+C’) 


ge b sin A’ 
sin B 


The area is computed with the formula 


1 


Area = = ac sin B 


This program works in any angular mode. However, if in degrees, decimal 
degrees are assumed. 
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OUTPU 


a 














Solve the triangle 














If B<90° and b <c, find al- 





ternate solution 

















* The stack must be maintained 








at these positions. 


Example: 
Given the following two sides and non-included angle: 


E B= 423° 


b=25.6 
c= 32.8 


Solve the triangle. 


Solution: 


Since B is less than 90° and b < c, two sets of solutions exist. 


C =59.58° 
A=78.12° 
a= 37.22 
Area = 410.85 
C’ = 120.42° 
A’ = 17.28° 
a’ = 11.30 


Area’ = 124.68 
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TRIANGLE SOLUTION a, b, c 


Given three sides of a triangle this program solves the triangle for the 
remaining parameters by the following formulas: 


2452. ¢2 
coors (Pat) 
cos rab 


p= sin-t (P88) A=sint (280€) 
c c 


This program also computes the area by the following formula: 
Area = +/s(s - a) (s ~ b) (s —c) 
1 
where s = — (atb+c) 


Reletter if necessary to make c the largest side. The program works in any 
angular mode. However, if in degree mode decimal degrees are assumed. 
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REGISTERS 








TOO 
RCL 1 
RCLO 
RCL 2 
RCLO 
RCL 3 


xs 
+ <|” 


N 
a 
oO 
oc} + 
N 
So 
st 
N 





2402 | RCL2 
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INPUT 
srer| INSTRUCTIONS DATA/UNITS 
Key in program 
Store sides (c is the largest) 


Solve the triangle 


























If only the area is needed: 

















* The stack must be maintained 





at these points. 


Example: 
Let a= 5.43, b = 10.46, c = 14.87 


Solution: 

C = 136.37° 
B= 29.04° 
A= 14.59° 
Area = 19.60 
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TRIANGLE SOLUTION a, A, C 


A c 


b 


Given two angles and an opposite side this program solves the triangle for the 
remaining parameters by the following formulas: 
B=2sin"' 1 -(A+C)=7 radians -(A + C) = 180° -(A+C) 
= 200 grads - (A + C) 
_ asin B 
sin A 
— asind 
sin A 


The area is computed with the following formula: 


Area = > ab sin C 


The program works in any angular mode. However, if in degree mode all 
angles are assumed to be in decimal degrees. 
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INPUT OUTPUT 


Key in program 

Store a, A, and C a 
A 
c 

Solve the triangle 

* The stack must be maintained 

at these points. 
































Let a = 19.6, A= 40.25°, C = 61.06° 


Solution: 

B= 78.69° 

b = 29.75 

c= 2655 
Area = 255.11 
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TRIANGLE SOLUTION a, b, C 


A 5 c 


Given two sides and their included angle this program solves the triangle for 
the remaining parameters by the following formulas: 


c=Va? +b? -2ab cos C A=sin"! (2%) 
c 


B=2sin"' 1 -(A+C) =m radians - (A + C) = 180° -(A+C) 
= 200 grads -(A + C) 


The area is calculated by 


Area = > ab sin C 


Reletter if necessary, to make a the smaller of a and b. 


This program works in any angular mode. However, if in degrees decimal 
degrees are assumed. 
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RCL 2 
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RCL 2 
RCL 3 
vx 
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RCL 1 
RCL 3 
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| 1405 | 
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* The stack must be maintained ms ses ee =e a 
|| [ 


Let a = 146, b = 227, C = 31.49° 


Solution: 

c= 127.76 

A= 36.65° 
B=111.86° 
Area = 8655.86 
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TRIANGLE SOLUTION a, B, C 


A 5 c 


Given two angles and their included side this program solves the triangle for 
the remaining parameters by the following formulas: 


A=2sin™' 1 -(B + C)=7 radians - (B + C) = 180° -(B+C) 
= 200 grads - (B+ C) 
_ asinB 
sin A 
_ asin C 
sin A 


The area is found using the formula: 


a? sin BsinC 


Me sin BEC) 


The program works in any angular mode. However, if in degrees the program 
assumes decimal degrees. 
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* The stack must be maintained 
at these points. 





Example: 
Let a = 20.96, B = 64°32’, C = 35°06’. 


Solution: 

First convert B and C to decimal degrees. 
A = 80.37° 

b=19.19 

c= 12.22 

Area = 115.66 
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HYPERBOLIC FUNCTIONS 


This program evaluates the six hyperbolic functions by the following 


formulas: 
xX =X 
. e*-e 
1. sinh x= 
2 


Xx ~X 
(Ce oh <z 
cosh x = —-_———. 


x ~xX 
ev =e 

tanks = —-——— 

e” re 


csch x = a (x #0) 
sinh x 


. sech x= 
cosh x 


coth x = ———_ 
tanh x 


gZ “n 
wz 


zs 
=-|> (rol + 
xX oO] |R ° x Z Pas 
° o kA a< 
a S Z 
cH 


DISPLAY 


REGISTERS 


xz 


F 00 | 
| 1 | 
oe | 
os | 
04 | 
05 | 
06 | 
Fos | 


7 
| 15:22 | 
mae 
| 13.00 | 
sina 
| 15:22 | 
| 51 | 
| 13.05 | 
aes 3] 
| 15-22 | 
aa 
Poet 
mee 
| 1473 | 


~s 
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14 73 
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LASTx 
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13 


rc 
> a nw = 
N N z 
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Examples: 
1. sinh 2.5 = 6.05 


cosh 3.2 = 12.29 
tanh 1.9 = 0.96 


sech (~.25) = 0.97 


2 
3 
4. csch 4.6 = 0.02 
5 
6 


coth (-2.01) = -1.04 
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INVERSE HYPERBOLIC FUNCTIONS 


This program evaluates the inverse hyperbolic functions by the following 


formulas: 


sinh”! x = In [x + (x? + 1)?] 
2. cosh™! x=In [x +(x? - 1)*] 


1. 


x21 





0<x<l 


a - 
ms] Oo ~~ 
> 
= a 
| 
a 
n 


5 8|8 8 
: g = =/2|" 
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Ce eo joe 


Ces DD Ld | 





Example: 


1, 


2 
3 
4. 
5 
6 


sinh™ (2.4) = 1.61 
cosh™! (90) = 5.19 
tanh~! (-.65) = -0.78 
esch™! (2) = 0.48 
sech™! (.4) = 1.57 
coth™! (3.4) = 0.30 





Photograph courtesy of NASA. 





HP-25 Program Form 


Ta a 
Switch to PRGM mode, press [f] [pram] , then key in the program. 
































newLerr Bl packarD 



























































INDEX 


Accumulated interest/remaining 
balance 32 

Area by double meridian distance 134 

Area of triangles 143-157 


Base conversions 22-25 


Calendar 49 
Chi-square evaluation 118 
Combinations 114 
Complex arithmetic 15 
Complex functions 18 
Compound amount 41 
Coordinate translation and rotation 138 
Course planning 61 
Covariance and correlation 
coefficient 101 
Cross product, vector 26 
Curve fitting 87-100 


Day of week, days between dates 49 

Determinant and inverse of 2 x 2 
matrix 20 

Discounted cash flow 46 

Distributions 105-109 

Dot product, vector 28 


Exponential curve fit 92 


Factorials 110 
Field angie traverse 129 


Great circle, plotting 62 
navigation 72 


Hyperbolics 158-161 


Interest rate, mortgage loan 39 
Internal rate of return 46 
Inverse, from coordinates 136 
hyperbolics 160 
normal integral 108 
ofa2x2 matrix 20 


Linear interpolation 85 
regression 87 
Logarithmic curve fit 95 


Moments and skewness 103 
Moon landing simulator 52 
Mortgage loan 32-40 


Net present value 46 
Newton’s method 76 
Nimb 55 
Normal distribution 105 
Numerical integration 81 
solution to differential equation 83 


One sample test statistics for the 
mean 127 


Periodic savings 44 
Permutations 112 
Plotting/graphing 7 
Power curve fit 98 
Probability 110-117 


Quadratic equation 12 


Random numbers 116 
Rhumbline navigation 65 


Sight reduction table 70 
Simultaneous equations 30 


t statistic, for two means 124 
paired 121 

Teach arithmetic 57 

Test statistics 118-128 

Traverse, field angle 129 

Triangle solutions 143-157 


Vector operations 26-29 
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